Stability of edge states in strained graphene 
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Spatially inhomogeneous strains in graphen e ca n simulate the effects of valley-dependent magnetic 
fields. As demonstrated in recent experiments^! the realizable magnetic fields are large enough to 
give rise to well-defined flat pseudo-Landau levels, potentially having counter-propagating edge 
modes. In the present work we address the conditions under which such edge modes are visible. We 
find that, whereas armchair edges do not support counter-propagating edge modes, zigzag edges do 
so, through a novel selective-hybridization mechanism. We then discuss effects of interactions on 
the stability of counter-propagating edge modes, and find that, for the experimentally relevant case 
of Coulomb interactions, interactions typically decrease the stability of the edge modes. Finally, we 
generalize our analysis to address the case of spontaneous valley polarization, which is expected to 
occur in charge-neutral strained graphene 1 ^'. 



I. INTRODUCTION 

Among condensed-matter systems, graphene is unique 
in being a flexible two-dimensional membrane whose 
electronic properties are tunable through deformations 
and strains. Spatially varying deformations affect the 
band structure of graphene by modulating the hopping 
amplitudes between lattice sites. Indeed, certain spa- 
tially varying deformation patterns can mimic the ef- 
fects of uniform "pseudo-magnetic fields," which have 
opposite signs in the two low-energy valleys^ (i.e., in 
the neighborhoods of the K and K' Dirac points) of 
graphene. (Unlike real magnetic fields, therefore, pseudo- 
magnetic fields preserve overall time-reversal invariance.) 
To date, such pseudo-magnetic fields have been realized 
using two distinct experimental approaches^! in both 
experiments, the fields realized were strong enough to 
drive the electronic structure deep into the quantum 
Hall regime in each valley, and — in line with theoret- 
ical predictions — the electronic structure measured by 
scanning-tunneling microscopy (STM) was clearly seen 
to consist of well-spaced pseudo-Landau levels (PLLs) 
in each valley. While the PLLs in each valley, consid- 
ered separately, have nontrivial topological invariant 
these are opposite for the two PLLs in the K and K' 
valleys; thus, the fate of various topological features in 
the full system is not completely understood. 

In the present work, we address the nature of one of the 
crucial topological features of strained graphene, namely 
its edge states. We first show, by analyzing the non- 
interacting system, that the Landau levels of strained 
graphene are very sensitive to the distinction between 
zigzag and armchair edges. For zigzag edges, the Landau 
levels (though flat in an infinite system) acquire an ap- 
preciable linear dispersion due to their hybridization with 
the pre-existing, non-topological "surface" states. (We 
shall use this term to distinguish these states from the 
topological "edge" states.) On the other hand, for arm- 
chair edges, where such boundary states do not exist, the 
Landau levels do not disperse even near the edge; this is 
in sharp contrast with the case of a real magnetic fielcP^ 



as we discuss, the physics behind this difference is that 
the lowest PLLs in the two valleys of strained graphene 
live on the same sublattice, whereas in the presence of a 
real magnetic field the corresponding Landau levels live 
on opposite sublattices. Thus, on a typical sample with 
rough edges, the edge states are expected to be localized 
in the zigzag regions and as a consequence should not 
contribute to transport; however, two-terminal measure- 
ments on zigzag nanoribbons should reveal the presence 
of multiple edge states. We then address the possibil- 
ity that Luttinger-liquid effects stabilize the edge states 
against disorder, but find that they do not for the case of 
a Coulomb interaction. Finally, we generalize our anal- 
ysis to the case of edge states in spontaneously valley- 
polarized states, and argue that, in general, these should 
not exhibit protected edge states. (One further mecha- 
nism for the destabilization of edge states, for instance 
in the experiments of Ref. [TJ is their hybridization with 
the Dirac sea in the surrounding, non-strained region; 
however, for undoped or weakly doped graphene this hy- 
bridization is presumably weak, owing to the vanishing 
density of states at the Dirac point.) 

Our paper is organized as follows. In Sec. [IT] we intro- 
duce the microscopic Hamiltonian for strained graphene 
that was used in our numerical work, as well as an effec- 
tive low-energy description that enables us to arrive at an 
analytic understanding of the effects discussed here. We 
then consider the physics of an armchair edge in Sec. Ill 
and that of a zigzag edge in Sec. [TVj In Sec. [V] we extend 
our analysis to include interaction effects and their influ- 
ence on edge stability. In Sec. |VI| we generalize our ar- 
guments, in the specific case of an armchair edge, to the 
case of a valley-polarized state. Finally, Sec. |VII| sum- 
marizes our results and discusses possible experimental 
signatures of edge-state physics in strained graphene. 



II. MODEL AND BULK PROPERTIES 

In what follows, we shall chiefly consider the fol- 
lowing noninteracting, nearest-neighbor Hamiltonian for 
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Strain generates a pseudo vector potential given by 
A° x + iA y = J2 a =i 2 3 St a (r)e ±iK - & a near the Dirac points 

±KP1. Note that E =i, 2 ,3 <^a(r)e ±,;K S ° is complex be- 
cause the nearest-neighbor hoppings are not symmetric 
under inversion. The real part of the strain gauge field A® 
is the same in both valleys and so couples to Qo = CTot 
and can be gauged away assuming time-reversal symme- 
try holds; whereas the imaginary part iAy, has opposite 
sign in the two valleys and couples witlP^l Q 3 = <j t 3 
leading to the valley-dependent magnetic fields realized 
in the experiments of Ref. [TJ 

For the purposes of the present work, we shall consider 
ribbon geometries, in which the strain-induced field is 
taken to realize the Landau "gauge" (A = (0, Bx)). A 
concrete lattice realization of this gauge field is shown in 
Fig. [I] In this realization of strain, exploiting translation 
invariance in y direction, the Dirac Hamiltonian can be 
written in the form: 



FIG. 1. Schematic pattern of change in hopping matrix el- 
ement on red bonds and varying along x direction as 8ts = 
evFBx which leads to the pseudo- vector potential A x = 
and A y — Bx. The dotted green square is the corresponding 
four cite unite cell and the green arrows are primitive lattice 
vectors. 



strained graphene: 
Ho = J2 12 (< + ^a(r l ))(a t (r t )6(r J + ^ a ) + ^. C .), (1) 

r 4 a=l,2,3 

where St a (ji) is, the strain-induced, nearest neighbour 
hopping amplitude modulation between the j4-sublattice 
site at Yi and the -B-sublattice site at + S a of the bi- 
partite honeycomb latticeP. The bond 8 a connects any 
A-sublattice atom to its three nearest neighbors on the 
B-sublattice. In the absence of strain, the low-energy ex- 
citations have a linear dispersion around the two Dirac 
points at momenta ±K with K = (47r/3v3cio) e x, o-o be- 
ing the carbon-carbon bond lengtrP. Near the Dirac 
points K and K' one can write the Bloch states in 
terms of a four-component spinor, as follows: \P = 
C0a,K) ipB,~K, V'AK'j ipB,K') where the first index denotes 
the component of the wavefunction on the A(B) sub- 
lattice of the honeycomb unit cell, and the second in- 
dex denotes the component of the state that is associ- 
ated with the K (K') valley. (In what follows, we shall 

write Va,k(K') = Wr(K') and ^b,k(K') = %(K')> in or- 
der to make contact with the standard notation for Dirac 
fermions.) The low energy effective Hamiltonian close to 
the Dirac points then reads as: 



Ho = v F 



-id x r 3 ai +p y T a 2 ~ -Bxt 3 <j 2 



(3) 



Ho = v F [p x T x + PyTy] 



(2) 



where T x — r 3 a\, Y y = to<?2, v f is the Fermi velocity, 
and the a and r operators are Pauli matrices acting on 
sublattice and valley indices respectively. 



where v F is the Fermi velocity, and the a and r oper- 
ators are Pauli matrices acting on sublattice (u and v) 
and valley (K and K') indices respectively. As this strain 
pattern respects translation invariance in the y direction, 
it is suitable for a strained ribbon extended in the y direc- 
tion. Notice that the Hamiltonian for valley K associated 
with p y is the same as the Hamiltonian for valley K' with 
momentum — p y . 

In this Landau "gauge", the wavefunctions in the 
zeroth PLL have the four-component form ^^py = 

(00,^,0,0,0) and ^ y = (0,0,^^,0), respectively, 
where 4>o,p y is the mth Landau orbital in the lowest (non- 
relativistic) Landau level, namely <f> Py cx exp(— ip y y— {x — 
Py^li) 2 /2^m) where Im is the magnetic length associated 
with the strain induced pseudo-magnetic field. In the 
higher PLLs, the wavefunctions take the form ^^ Py = 

(<f>n,py,<f>n-i,p v ,Q,0) and * K Py = (0,0,0* ,0*_ liPv ). 
These forms should be contrasted with those for graphene 
in a real magnetic field; in the case of a real field, the 
Landau level wavefunctions in the two valleys have oppo- 
site sublattice structure, whereas in the case of a strain- 
induced field, the wavefunctions have the same sublattice 
structure. 

In our treatment of the nonintcracting problem we 
have suppressed the physical spin, and thus ignored the 
(weak) intrinsic spin-orbit coupling in graphene. How- 
ever, the spin index can be included trivially. We shall 
consider the consequences of the physical spin in Sec. |VJ 
as interaction effects in the spinless and spinful cases are 
different; moreover, a combination of interaction effects 
and spin-orbit coupling was argued to lead to stabiliza- 
tion of fractional phases^. (The spinless situation can be 
experimentally realized by applying a large field paral- 
lel to the graphene sheet, and thus spin-polarizing the 
electrons.) 
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III. ARMCHAIR EDGE 

Having discussed the bulk properties, we now turn to 
the case of a graphene ribbon with armchair edges along 
the y direction. Owing to the translational invariancc 
along y as well as our choice of the Landau "gauge" 
the momentum in y direction is a good quantum num- 
ber. On the other hand, in the ribbon geometry, trans- 
lation invariance in x direction is broken. For our ana- 
lytic calculations we consider the case of a semi-infinite 
graphene ribbon covering the region x < 0, with an arm- 
chair boundary at x = 0. 




If we make the transformations uk' -uk' and 
«K' — > —vjt>, the Hamiltonian in the two valleys reads 



as: 



H K , = v F 



-id x ai + p y a 2 Bxa 2 

c 

-id x ai ~p y <J 2 ~ -Bxa 2 



(4) 
(5) 



In order to get the spectrum, we square these Hamilto- 
nians, leading to: 



-dl 

-dl- 



Py 



Py 



Bx - ~Bao 



-Bx 



-Ba^ 



(6) 
(7) 



which is defined for x < 0. In the zeroth Landau level, in 
which (for a strain-induced field) the wave functions in 
both valleys are based on the same sublattice, so that the 
a matrix is trivial. The boundary condition then reduces 
to the wave function vanishing on the last row. In terms 
of the sublattice wave functions this reads as: 



-1 -0.5 0.5 1 

Py 



(a) 




-1 -0.5 0.5 1 

Py 



(b) 

FIG. 2. Energy dispersion for a graphene nanoribbon with 
armchair edges as a function of the momentum along the edges 
(p y ); (a) in the presence of a strain- induced pseudo magnetic 
field and (b) in the presence of a real magnetic field. The 
center of mass position of the Landau orbitals is correlated 
with their momentum p y i.e. p y — is localized in the center 
of the ribbon and as \p y \ increases, they get closer to the edge. 
For the strain induced magnetic field, there is no dispersion 
in the zeroth PLL, whereas for real magnetic field there are 
dispersing edge states associated to the topological character 
of the bulk (i.e., the quantum Hall effect). 



uk(p v ) = UK< (Py) (8) 

If we change x — > —x for K', the two equations ^ and 
([7]) become identical. We can join the two equations and 
the boundary condition above will then correspond to 
continuity of the wave function at x = 0. The Hamil- 
tonian then reduces to the simple harmonic oscillator 
Hamiltonian and the boundary question is trivially sat- 
isfied. The salient property of this Hamiltonian is the 
fact that the energy eigenvalues are independent of p y , 
and thus the states near the edge do not disperse. This 
is i n sh arp contrast with the case of a regular magnetic 
fielcP^, in which states near the edge do disperse. 

One can understand the nondispersing nature of the 
edge states intuitively as follows, for the case of a semi- 
infinite ribbon with an edge at x = 0. Because both the 
(K,k y ) and (K',k y ) wavefunctions are located on the 
same sublattice, and their guiding centers are at k y and 
— k y respectively. One can always construct equal- weight 
superpositions of the two wavefunctions, which automat- 
ically satisfy the boundary condition. (By contrast, in 
the case of a regular magnetic field, the two states live on 
opposite sublattices and therefore cannot be superposed 
to meet the boundary conditions; instead, the bound- 
ary condition must be satisfied by introducing a slowly 
varying envelope functional, and this costs additional ki- 
netic energy for states near the edge.) Note that these 
arguments extend straightforwardly to the case of higher 
PLLs. 

As Fig. [2ja) shows, a numerical computation of the 
band structure is consistent with the argument above. 
An important implication of this argument is that there 
should be no counter-propagating edge modes on arm- 
chair edges; i.e., if the chemical potential lies between 
two PLLs, it should not cross any states at all. 
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IV. ZIGZAG EDGE 



We now turn to the case of nanoribbons with zigzag 
edges. Even in the absence of a magnetic field, these 
edges support non-dispersing modes localized at the 
edges of the sample; these "surface" states occupy the 
A sublattice at one end of the system, and the B sublat- 
tice at the other end. Because these states are confined 
on scales smaller than the magnetic length, they are ex- 
pected to be unaffected by a magnetic or strain fielcP^. 
In the case of a real magnetic field, the dispersion of 
the topological edge states^ occurs chiefly by intrinsic, 
universal means, and is therefore only slightly modified 
by the presence of zigzag surface states. The case of a 
strain-induced magnetic field is fundamentally different 
in this respect: as we have argued above, there is no in- 
trinsic dispersion mechanism, and thus the primary con- 
tribution to the Landau level dispersion is the differential 
hybridization between Landau levels and the zigzag edge 
(or "surface") states. 

We now address the effects of this hybridization, fo- 
cusing at first on the lowest PLL. In the lowest PLL, 
all Landau orbitals are situated on either the A or the 
B sublattice. Thus, they hybridize only with the zigzag 
surface states at one edge of the system; the strength of 
this hybridization depends on the distance of the guid- 
ing center of the Landau orbital from the "Hybridizing" 
edge. 

Assuming the hybridizing edge is located at x = 0, 
this hybridization changes the energy of each Landau or- 
bital by a factor ~ exp[— (jp v Im) 2 }- Because this energy 
shift is larger for levels situated near the edge, it causes 
these states to disperse. Moreover, because of the valley- 
dependence of the pseudo-magnetic field, the states in the 
K valley that are near the edge have opposite momenta 
to the corresponding states in the K' valley. Near the 
"non-hybridized" edge, by contrast, there is no hybridiza- 
tion and thus no dispersion at all. The size of this effect 
is quite strong in mesoscopic systems (such as all exist- 
ing experimental realizations of strain-dependent gauge 
fields), as one can see from diagonalization on graphene 
nanoribbons (Fig. |3(aj| ). 

Two important observations can be made about this 
surface-induced dispersion effect. 

(1) The hybridization with the zigzag edge state plays 
a role similar to an electric field perpendicular to the 
edges of the ribbon; thus, the dispersion generated by the 
zigzag surface states can presumably be either enhanced 
or compensated by the application of such an electric 
field. In particular, the number of surface states cross- 
ing the chemical potential might be detectable by these 
means. (2) Unlike armchair edges, zigzag edges with sur- 
face states do not mix the valleys. This can easily be seen 
from Fig. 3(a) — the level crossings between PLLs from 




(a) 




opposite valleys are not avoided, as they would be in the 
presence of valley mixing. The physical reason for this is 
that the Landau orbitals in the two valleys — since they 
are counter-propagating — couple to orthogonal, counter- 



00 



FIG. 3. Energy dispersion for a graphene nanoribbon with 
zigzag edges as a function of the momentum along the edges 
(p y ); (a) in the presence of inversion symmetry and (b) in the 
presence of a inversion symmetry breaking perturbation. 



propagating linear combinations of the zigzag edge states. 
Thus, for a clean zigzag edge, there should be counter- 
propagating edge states, which can be detected (e.g.) in 
two-terminal transport measurements. Note that these 
counterpropagating edge states can be gapped out by 
adding a valley-mixing term such as the Kekule distor- 
tion (see Fig. |3(b)| ). 

In practice, perhaps the most promising candidate for 
realizing a system with a clean zigzag edge is molecu- 
lar graphene. In this system, one can tunably create a 
Kekule distortion 2 and thus study the transition between 
Fig. [3(a)] and Fig. [3(b)] 
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V. STABILITY UNDER INTERACTIONS 



The considerations discussed so far apply to the case 
of noninteracting edge modes. We now address the is- 
sue of whether the edge modes continue to carry current 
in a disordered interacting system, using the approach 
of Ref. 1 -^. (Similar questions were addressed for quan- 
tum spin-Hall insulators with multiple edge modes in 
Ref. [T^l) It is clear that, in the absence of interactions, 
backscattering due to disorder localizes the edge modes, 
leading to an insulating edge. However, there are certain 
regimeifSE^ in which interactions prevent the localiza- 
tion of edge states, by making backscattering processes 
"irrelevant" in the renormalization-group sense. (Quali- 
tatively, this effect can be understood in terms of inter- 
actions screening out the impurity potential; this would 
happen, e.g., if the electron-electron interactions were at- 
tractive™.) 

For repulsive interactions, as discussed in Ref. [ITJ the 
low-temperature behavior of the conductance is given 
by the strength of the density-density interaction be- 
tween the left-moving and right-moving edge states (i.e., 
forward-scattering) compared with back-scattering pro- 
cesses involving momentum-exchange. For the conduc- 
tance to stay finite as T — > 0, the above criterion implies 
that the Fourier components of the interaction potential 
satisfy V(0) < |V(2|G|), where G is a reciprocal lattice 
vector. (This is an approximate expression; the precise 
momentum transfer involved in backscattering is the mo- 
mentum difference between the two states on the edge, 
which is ~ |G| + 2k p, which we approximate as |G| us- 
ing the assumption that the chemical potential is near the 
Dirac point.) This criterion is not met by realistic inter- 
action potentials. However, it has been showrpJJ that for 
V(0) < 2V(2|G|) (as in, e.g., a contact potential, or a 
sufficiently short-ranged potential achieved via gating), 
the conductance exhibits non-monotonic behavior and 
increases with decreasing temperature until a crossover 
scale 



_Dexp 



2F(2|G|)-y(0) 
>(2|G|)(2V(0)-V(2|G|)) 



(9) 



where D is the bandwidth of the Dirac band (for molec- 
ular graphene this is ~ 200 meV). The physics of this 
crossover can be explored in strained graphene by tuning 
the range of the inter-electron interaction via gating. 



VI. EFFECTS OF VALLEY POLARIZATION 

Thus far, we have considered the edge states of strained 
graphene under the assumption that time-reversal invari- 
ance is preserved. We now turn to situations in which 
time- reversal invariance is spontaneously broken; such 
situations naturally arise in the case of half-filled PLLs 3 , 
in which the valley-polarized state naively appears to be 




CD W order 
at domain wall 



FIG. 4. Corbino disk with broken valley degeneracy 



a topologically protected state. We argue in this section 
that such states do not generically have topologically pro- 
tected edge modes. 

Within the low-energy theory, there are two valley- 
polarizing perturbations, given by t z oq (i.e., simple 
valley-polarization) and t z o z (i.e., the Haldane mass 
terrrP^). The former opens up a gap in any half-filled 
PLL, whereas the latter opens up a gap only in the zeroth 
PLL. We shall consider these terms separately, beginning 
with the simple valley-polarization term. As discussed in 
Ref. 1 101 this term anticommutes with two charge-density- 
wave perturbations that take the form t x (Jq and t v o-q in 
the low-energy theory. Being spatial modulations, these 
terms are naturally generated at sharp edges; thus, the 
valley polarization gap in the bulk can be rotated into the 
valley-mixing gap at the edge without closing any gaps, 
thus leading to the absence of edge states. We only con- 
sider the case of uniform charge-density-wave terms on 
the edge, and neglect the question of whether disordered 
or local charge-density-wave terms can gap out the edges. 
One expects this mechanism to prevent edge states from 
being present at half-filling of any PLL, as long as (a) the 
edges mix the valleys, or (b) the charge-density-wave is 
dynamically generated by interactions at the edge. It is 
of course possible that very smooth edges, which do not 
mix the valleys, can support chiral edge states, but the 
fragility toward valley-mixing indicates that these edge 
states are not in fact protected. 

Before we proceed to the case of the zeroth PLL, we 
briefly explain how our remarks above can be related 
to the conventional flux-insertion argument!^, which was 
usecP to argue for topologically protected edge modes in 
half-filled PLLs. Consider half-filled PLLs in a Corbino 
disk geometry as shown in Fig. [4j in which the left half of 
the disk has a perturbation of the form r 2 and the right 
half has a perturbation of the form — r z ; this situation, 
which corresponds to magnetic domains, should typically 
arise in experiments, given the Ising nature of the valley 
polarizatiorP. We assume that the gauge realized by the 
strain has the property that p y = is near the center of 
the disk, and then add a perturbation of the form t x in 
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the medial region of the Corbino disk. (This assumption 
can easily be relaxed, if one then chooses a perturbation 
that is modulated along y.) Now suppose one inserts a 
(real) magnetic flux through the center of the Corbino 
disk: if the regions having a t z perturbation were indeed 
topological, one would expect the flux insertion to pump 
two charges from the edges of the disk to its center (given 
the opposite nature of the strain-induced magnetic field 
in the two regions). However, this is impossible because 
the center of the disk is fully gapped and therefore in- 
compressible; thus, we conclude that the flux insertion 
argument fails, and that half-filled PLLs do not in gen- 
eral have a quantized Hall conductance. 

The argument given above holds for all the PLLs; how- 
ever, an additional subtlety arises in the case of the ze- 
roth PLL, owing to the possibility of a Haldane mass 
gap, T z afi**. In any PLL other than the zeroth, this term 
does not open a gap. However, within the zeroth PLL, 
it does open up a gap, which leads to the following ap- 
parent paradox. As the Haldane mass term generates 
a gap both with and without the strain-induced field, 
one can imagine adding it before adding strain; in this 
case, the gap generated is a topologically nontrivial gap, 
accompanied by chiral edge states. The strain-induced 
field does not compete with this gap (because the lowest 
PLL is gapped out by the Haldane mass), thus this topo- 
logically trivial gap and the corresponding edge states 
must continue to exist even in strained graphene with 
a Haldane gap. On the other hand, the Haldane gap — 
when projected onto the zeroth PLL — can be rotated into 
the charge-density- wave gaps discussed in Ref. [TUJ thus, 
its projection is topologically trivial and cannot generate 
edge states. Our numerical diagonalization studies, plot- 
ted in Fig.[5j suggest that the resolution is as follows — the 
topologically protected Haldane edge state is connected 
to higher-order PLLs, and crosses the zeroth PLL with- 
out hybridizing with it. Thus, although these edge states 
give rise to a robust, quantized Hall effect, the mecha- 
nism is unrelated to the low-energy PLL structure. In- 
deed, in the limit where the Haldane mass term is much 
smaller than the other perturbations, we find that these 
edge states move far away from the Dirac points, demon- 
strating that they are unrelated to the PLLs. In ad- 
dition, when the electron density is tuned to be at the 
charge neutral point, and the valley polarization (hence 
the Haldane mass) is generated spontaneously by inter- 
actions, the interaction effects are operative nominally 
only for the half-filled lowest PLL, but not for the other 
PLLs that are totally empty or occupied. I.e., the Hal- 
dane mass is selectively generated for the lowest PLL, 
but not for the other PLLs. We would then expect there 
is no edge state at all in this case. 



VII. CONCLUSIONS AND OUTLOOK 

In the present work we have addressed the properties 
of the edge states of strained graphene in the quantum 
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FIG. 5. Pseudo-Landau level (PLL) structure in the presence 
of a large Haldane mass gap. As one can see, the topologically 
protected edge modes guaranteed by the Haldane mass term 
do not mix with the PLLs, but appear to originate at much 
higher energies. 



Hall regime. We have argued that the edge physics is not 
universal — owing to the lack of topological protection — 
but is strongly dependent on the nature of the edge. 
Whereas armchair edges do not support edge states 
at all, zigzag edges are expected to support counter- 
propagating edge states. Furthermore, in the case of 
zigzag edges, the hybridization between the PLL states 
and the zigzag "surface" states gives rise to dispersion 
of the PLLs; for short-range interactions, these counter- 
propagating edges are expected to manifest themselves 
via nonmonotonic temperature-dependence of the con- 
ductance. Finally, we considered the case of valley- 
polarized edges (e.g., a quantum Hall ferromagnet in 
charge-neutral strained graphene) and argued that, in 
general, valley-polarization does not imply a finite Hall 
conductance, because (within the zeroth PLL) valley po- 
larization can be continuously transformed into a charge- 
density wave, which is evidently non-topological. 

Our results have several experimental implications. 
Most notably, the difference between the dispersion near 
zigzag and armchair edges can easily be detected using 
the spectroscopic methods of current experiments. In ad- 
dition, transport experiments would — assuming the con- 
tribution of the PLLs could be isolated — provide several 
clear signatures. For example, for strained graphene rib- 
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bons having clean zigzag or armchair edges, the strik- 
ing difference between the former case (with multiple 
countcrpropagating edge modes) and the latter (no edge 
modes at all) should be easily detectable via two-terminal 
measurements. Similarly, a nonmonotonic temperature- 
dependence of the two-terminal conductance would pro- 
vide a signature of interaction effects near the edges. Fi- 
nally, either spectroscopy or transport can address our 
prediction that, for half-filled PLLs, the regions between 
magnetic domains should (a) contain no propagating 
modes, and (b) exhibit spontaneously broken transla- 



tional symmetry. 
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